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Abstract 

We investigate the evolution of a scalar field propagating in Reissner-Nordstrom 
Anti-de Sitter spacetime. Due to the characteristic of spacetime geometry, the radiative 
tails associated with a massless scalar field propagation have an oscillatory exponential 
decay. The object-picture of the quasinormal ringing has also been obtained. For 
small charges, the approach to thermal equilibrium is faster for larger charges. How- 
ever, after the black hole charge reaches a critical value, we get the opposite behavior 
for the imaginary frequencies of the quasinormal modes. Some possible explanations 
concerning the wiggle of the imaginary frequencies have been given. The picture of the 
quasinormal modes depending on the multipole index has also been illustrated. 
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1 Introduction 



The study of wave dynamics outside black holes has been an intriguing subject for the 
last few decades (for a review, see Q). In virtue of previous works, we now have the 
schematic picture regarding the dynamics of waves outside a spherical collapsing object. 
A static observer outside the black hole can indicate three successive stages of the wave 
evolution. First the exact shape of the wave front depends on the initial pulse. This stage 
is followed by a quasi-normal ringing, which carries information about the structure of the 
background spacetime and is believed to be a unique fingerprint to directly identify the black 
hole existence. Detection of these quasinormal modes is expected to be realized through 
gravitational wave observation in the near future [|]]. Finally, at late times, quasinormal 
oscillations are swamped by the relaxation process. This relaxation is the requirement of 
the black hole no-hair theorem [Q. Besides the dynamical mechanism of shedding the 
perturbation hair near black hole event horizon is of direct interest to the problem of stability 
of Cauchy horizons [[|. 

The mechanism responsible for the relaxation of neutral external perturbations was first 
exhibited by Price ||. Studying the behavior of a massless scalar field propagating on a 
fixed Schwarzschild background, he showed that for a fixed position the field dies off with a 
power-law tail. The behavior of neutral perturbations along null infinity and future event 
horizon was further studied by Gundlach, Price and Pullin || and similar power-law tails 
have been obtained. These results were later confirmed using several different techniques, 
both analytic and numerical || |7|, and were generalized to Reissner-Nordstrom (RN) 
background |5], |§. The application of linear approaches is encouraged by numerical analysis 
of the fully nonlinear dynamics of the fields |T0| , [IT| , which indicates the same late time 
pattern of decay. 

Extending the basic scenario to study the asymptotic evolution of charged fields and 
self-interacting massive scalar field around a RN black hole, Hod and Piran |}I2| , [EJ found 
that although usual inverse power-law relations of these fields present at timelike infinity and 
null infinity, along the future black hole event horizon this power-law tail is accompanied 
by an oscillatory behavior. Recently, Brady et al. |14|, |1"5|| studied scalar wave dynamics in 
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non-asymptotically flat exteriors of Schwarzschild-de Sitter and RN de Sitter black holes. 
Contrary to the asymptotically flat geometries, no power-law tails were detected in these 
cases. Instead the waves were found to decay exponentially at late times. For 1 = 0, they 
found that the field does not decay, but settles down to a nonzero constant. Moreover, for 
a field strongly coupled to curvature, they obtained that the wave function oscillates with 
an exponentially decaying amplitude. These observations support the earlier argument by 
Ching et al. |16| that usual inverse power-law tails as seen in asymptotically flat black 
hole spacetimes are not a general feature of wave propagation in curved spacetime. Besides 
some relation to the perturbative field, the relaxation process reflects a characteristic of the 
background geometry. 

It is of interest to extend this study to Anti-de Sitter (AdS) spacetime. In addition to 
three major aspects that the evolution of test-field is associated with, including the no-hair 
theorem, the stability of Cauchy horizon and direct evidence of the existence of black hole 
provided by quasinormal ringing, the recent discovery of the Anti-de Sitter/ Conformal Field 
Theory (AdS/CFT) correspondence makes the investigation in AdS black hole background 
more appealing. The quasinormal frequencies of AdS black hole have direct interpretation in 
terms of the dual Conformal Field Theory (CFT). In terms of the AdS/CFT correspondence 
171, [TjL [HI], the black hole corresponds to an approximately thermal state in the field theory, 



and the decay of the test-field corresponds to the decay of the perturbation of the state. 
The first study of the scalar quasinormal modes in AdS space was performed by Horowitz 
and Hubeny on the background of Schwarzschild AdS black holes in four, five and seven 
dimensions. They claimed that for large black holes both the real and imaginary parts of 
the quasinormal frequencies scale linearly with the black hole temperature. The time scale 
for approaching the thermal equilibrium is detected by the imaginary part of the lowest 
quasinormal frequency and is proportional to the inverse of the black hole temperature. 
Considering that the RN AdS solution provides a better framework than the Schwarzschild 
AdS geometry and may contribute significantly to our understanding of space and time, we 
generalized the study made in our previous work [|22| . We found that the charge in RN AdS 
black hole showed a richer physics concerning quasinormal modes and further information 
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on AdS/CFT correspondence. The bigger the black hole charge is, the quicker for their 
approach to thermal equilibrium in CFT. However these studies focused much on frequencies 
of the quasinormal modes, the object-picture of the evolution of test-field around the AdS 
background is lacking. 

The intention of this paper is to analyse in detail the wave propagation of massless 
scalar field in RN AdS spacetime. We will show that the direct picture of the evolution 
presents us with a perfect agreement on quasinormal frequencies with those obtained by 
using approximation method suggested in ]20|]. Moreover, our work will show object-pictures 
of the behavior of quasinormal modes as a function of the charge Q and of the multipole 
order of the field /. We will also address some discussions on highly charged background 
case. In addition, the relaxation process at the event horizon will also be discussed in our 
direct picture. We found that the decay has the pattern of oscillatory exponential tail. This 
result supports Horowitz's claim that there are no power-law tails at late times in AdS space. 
Some physical explanation related to this result will be given. 

2 Equations and numerical methods 

The Reissner-Nordstrom black hole solution of Einstein's equations in free space with a 
negative cosmological constant A = —3/R 2 is given by 



ds 2 



hdt 2 + /rMr 2 + r 2 dn 2 



(1) 



with 



^ r + r+ Q 2 Q 2 r 2 



r R 2 r r + r r 2 R 2 



(2) 



The asymptotic form of this spacetime is AdS. The mass of the black hole is 




(3) 



The Hawking temperature is given by the expression 



and the potential by 



Q_ 



In the extreme case r + and Q satisfy the relation 

ri + i? 2 



(5) 



(6) 



For a non-extreme RN AdS black hole, the spacetime possesses two horizons, namely 
the black hole event horizon r + and Cauchy horizon r_. For the extreme case where (|6]) is 
satisfied, these two horizons degenerate. The function h has four zeros at r + ,r_ and r 1; r 2 , 
where ri,r 2 are two complex roots with no physical meaning. The relations between these 
four roots are 



ri + r 2 = -(r + + r_) 

r 1 r 2 = R 2 + r + r^ + r 2 , + r 2 _. 



(7) 



In terms of these quantities, h can be expressed as 

1 
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Introducing the surface gravity associated with r; by the relation 
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These quantities allow us to write 

1 



/T 1 



+ 



2n r+ (r — r + ) 2K r _(r — rJ) 2n ri (r — r x ) 2n r2 (r — r 2 ) 



(10) 
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Combining the last two terms in (10), we express the transformation between r and the 
"tortoise coordinate" r* = J h~ l dr in the form 



2k 



+ 



m(r — r + ) ln(r — r_) 

R 2 [{r 1 r 2 ) 2 — r + r_rir 2 ] 



[r% - r+(ri + r 2 ) + rir 2 ][r 2 - r_(r x + r 2 ) + r 1 r 2 ] 
i? 2 [r + r_(ri + r 2 ) — r+r\r 2 — r_rir 2 ] 



y r 2 — 



r(ri + r 2 ) + rir 2 
rdr 



(11) 



[r\ - r + (r 1 + r 2 ) + rir 2 ][r 2 - r_(r x + r 2 ) + rir 2 ] J r 2 - r(r 1 + r 2 ) + rir 2 

In terms of t and r* we introduce null coordinates u = t — r* and v = t + r* so that 
the future black hole horizon is located at u = oo. Since the quasinormal modes of AdS 
space are defined to be modes with only ingoing waves near the horizon, we will pay more 
attention on the wave dynamics near the event horizon. 

Let us consider a massless scalar field $ in the RN AdS spacetime, obeying the wave 
equation 

□ $ = (12) 

where □ = g a/3 'V a Vi3 is the d'Alembertian operator. If we decompose the scalar field ac- 
cording to 

<f> = J2-Mt,r)Y lm {6,(f>) (13) 
, r 

then each wave function ipiir) satisfies the equation 

d 2 ipi d 2 i)i 



where 



Vi 



1(1 



1(1 + 1] 



dt 2 



ldh 

r dr 
r+ -f 



<9r* 5 



(14) 



rl/R 2 + Q 2 /r + 2Q 2 



2 

B 2 



(15) 



The potential Vi has the same characteristic as that in Schwarzschild AdS black hole. 
It is positive and vanishes at the horizon, but diverges at r — > oo, which requires that $ 
vanishes at infinity. This is the boundary condition to be satisfied by the wave equation 
for the scalar field in AdS space. In terms of the radial coordinate r*, it is seen that when 
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r tends to infinity, r* tends to a finite constant (which we denote r* s ). It means that our 
region of interest in the (u — v) diagram is above the line v — u — 2r* s , as shown in figure 1. 
In this line (where r — > oo) we set $ = 0. 

The behavior of the potential differs quite a lot from that of asymptotically flat space 
and de Sitter space. As argued in [|1(|, it is this peculiarity that contributes to the special 
wave propagation as will be shown later. 

Using the null coordinates u and v, the equation (14) can be recast as 

- 4 g^M^ v) = Vtir^u, v) (16) 

in which r is determined by inverting the relation r*(r) = (v — u)/2. 

The two-dimensional wave equation (16) can be integrated numerically, using for example 
the finite difference method suggested in []5], ^TJ. Using Taylor's theorem, it is discretized as 

iPn = We + Ww - Ws - oudvVi I j hO(e) (17) 

where the points N, S, E and W form a null rectangle with relative positions as: 
N : (u + 5u,v + 5v),W : {u + 5u,v),E : (u, v + 5v) and S : (w, v). The parameter e is 
an overall grid scalar factor, so that 5u ~ 5v ~ e. 

Considering that the behavior of the wave function is not sensitive to the choice of initial 
data, we set Wi(u, v — v o) = and use a Gaussian pulse as an initial perturbation, centered 
on v c and with width a on u = uq as 



Wi{u = u Q ,v) = exp 



2a 2 



(18) 



The inversion of the relation r*(r) needed in the evaluation of the potential VJ(r) is the 
most tedious part in the computation. We overcome this difficulty by employing the method 



suggested in [H, 14 
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Figure 1: (left) Diagram of the numerical grid and the domain of interest. The black spots represent 
the grid points where the value of the field is known. The gray spots represent the grid points to be 
calculated. The points in the line u = u max are the results shown in this paper, (right) Detail of 
the previous diagram, showing the relative positions N, S, E and W. The gray point would be the 
next one to be calculated. 

After the integration is completed, the value ipi(u max ,v) are extracted, where u max is 
the maximum value of u on the numerical grid. Taking sufficiently large u max , ipi(u max ,v) 
represents a good approximation for the wave function at the event horizon (u — > oo), which 
carries information about the quasinormal modes for AdS space of our interest. It was 
observed in our numerical experiments that the plots obtained converged in a given range 
of v, as expected, but the rate of convergence varied with the parameters of the system. It 
should be noted that the term "convergence" in the last sentence refers to the limit of $ 
as u max is taken bigger and bigger. As commented, the region u max — > oo corresponds to 
the event horizon, in which we are interested. Some comments about the convergence of the 
numerical code, meaning the evolution of the wave function with the grid size, is made in 



appendix A. In order to compare our results here with those in J2(| ^2f , we fix R = 1 in the 
following. 



S 



3 Numerical results 



We now report on the results of our numerical simulations of evolving massless scalar field 
on a RN AdS black hole background. Taking Q — > 0, our results reflect the properties on 
Schwarzschild AdS background. 

3.1 Behavior of wave evolution for / = 

In the first series of numerical experiments, we choose the multipole index I = and examine 
the behavior of the test-field propagation with the increase of the charge of the background 
spacetime. 




-3.5 1.5 6.5 11.5 10 20 30 



v v 

Figure 2: Semi-log graph of the absolute value of the wave function for small values of v (left) and 
the late time behavior (right), with r + = 0.4 and Q = 0. 

The figure 2 displays the quasinormal ringing on the 4D Schwarzschild-AdS background. 
From plots using different values of r + we can read off values for the imaginary (ujj) and real 
(cor) parts of the frequency The results are shown in table 1. 
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Table 1: The lowest quasinormal mode frequency for four-dimensional Schwarzschild-AdS black 
hole for 1 = 0. The real part is uj r and the imaginary part is uj. 

The agreement with the frequencies calculated in is good. A small difference can be 
attributed to the different methods employed. 

We studied the late-time decay of the test-field on Schwarzschild-AdS black hole and the 
results are shown in figure 2. It was first predicted in j20j that the decay is always exponential 
and with no power-law tails. Through careful study, we got a result that supports their 
claim. We obtain that the late-time falloff is oscillatory exponential. Recall that for the 
AdS black hole, the potential diverges at infinity but vanishes exponentially near the black 
hole horizon. This boundary condition differs completely from that of asymptotically flat 
spacetime. Using Ching et al's argument fTE] , it is not difficult to understand the reason 
that usual power-law tail is replaced by exponentially decay. Since that potential also falls 
off exponentially on black hole horizon in de Sitter space, the pattern of decay exhibit here 
on AdS black hole event horizon is similar to that in de Sitter space [14], 0. However 
due to the waves bouncing off the divergent potential barrier at large r in AdS space, the 
oscillation appears in the exponential tail. With the late-time behavior in hand, it is possible 
to improve previous studies on Cauchy horizon stability problems in AdS black holes |23|, f24}j 
along the lines of Brady and Smith ||25|| . 
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Figure 3: Semi-log graphs of |^| wi/t r + = 0.4 and small values of Q. The extreme value for Q 2 
is 0.2368. 
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Table 2: The behavior of the lowest quasinormal mode frequency with the increase of the charge in 
RN AdS black hole background for I = 0. The real part is and the imaginary part is uj}. 

Figure 3 demonstrates the behaviors of the field with the increase of the charge in RN AdS 
black hole background. In table 2 we listed values of imaginary and real parts of quasinormal 
frequencies read from the plots. These frequencies together with the entire picture presented 
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in figure 3 agree perfectly with the result given in [p2| . Since the imaginary and real parts 
of the quasinormal frequencies relate to the damping time scale (ti = 1/uJi) and oscillation 
time scale (r 2 = l/u r ), respectively. We learned that as Q increase, uii increase as well, 
which corresponds to the decrease of the damping time scale. According to the AdS/CFT 
correspondence, this means that for bigger Q, it is quicker for the quasinormal ringing to 
settle down to thermal equilibrium. From figure 3 it is easy to see this property. Besides, 
table 2 and figure 3 also tell us that the bigger Q is, the lower frequencies of oscillation will 
be. This object-picture support the argument in that if we perturb a RN AdS black 
hole with high charge, the surrounding geometry will not "ring" as much and as long as that 
of the black hole with small Q. It is easy for the perturbation on the high charged AdS black 
hole background to return to thermal equilibrium. However this relation seems not to hold 
well when the charge is sufficiently big and near the extreme value described by relation (|6]). 
We will address this behavior later. 



3.2 Behavior of wave evolution with the increase of / 

We have so far discussed only of the lowest multipole index 1 = 0. In the following we 
show the wave dynamics on the RN AdS background for different multipole indices /. For 
simplicity we present here the results for Z = 0, 1, 2, but the basic characteristics of the 
dependence of the wave function with the multipole index have already been uncovered and 
are discussed in this section. 
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Figure 6: Semi-log graph of the wave function for Q 2 
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= 0.125 and r + = 0.4, with I = 0, 1, 2. 



Figures 4 to 6 exhibit a striking consistent picture with those given in [pOj , 22 1. Increasing 
/, the evolution of the test field experiences an increase of the damping time scale and 
a decrease of the oscillation time scale. These figures gave us an object-lesson on the 
evolution of the test field with the increase of multipole index. It is worth noting that the 
behavior shown here differs quite a lot from that of the asymptotically-flat case ||, where 
the perturbation settles down faster with the increase of /. These differences can be used 
to further support the argument that the last two processes of the wave evolution reflect 
directly of the background spacetime property. 



3.3 Quasinormal modes for highly charged AdS black hole 



In our previous work p2 |, we found that there is a numerical convergence problem when 
the charge Q approaches the extreme value. We claimed that the problem is related to the 
method we adopted there. Using the numerical strategy proposed in || to directly describe 
the wave dynamics, we can step further on this problem. 
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Figure 7: (left) Semi-log graph of the wave function with r + = 0.4, I = and several values of Q. 
The extreme value for Q 2 is 0.2368. (right) Detail of previous plot, with two values ofQ. 

The behavior we read from figure 7 is quite different from that exhibited in figure 3 and 
discussed in subsection 3.1. Here we see that over some critical value of Q, the damping time 
scale increases with the increase of Q, corresponding to the decrease of imaginary frequency. 
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This means that over some critical value of Q, the larger the black hole charge is, the slower 
for the outside perturbation to die out. This qualitative change in the characteristic of the 
imaginary frequency as Q increases agrees to the normal-mode frequencies of RN black hole 
described in [^f], |2?J . 

We found in the numerical calculation that the computation time becomes bigger when 
Q is large enough. This happens because we are forced to take a bigger u max , and also 
because the numerical routine which calculates the function r(r*) is more time consuming 
in this range of Q. Moreover, the precision cannot be improved for taking smaller grid scale 
factor when Q is close to the extreme value. Some plateaus appear in the late-time leading 
us to suspect that there may be a failure in the numerical method in this region. In spite 
of this, the convergence of the numerical code appears to be good even in this limit, as 
commented in appendix A. However the question whether there is a plausible explanation 
for the "wiggle" of the imaginary frequency as the charge is large enough and close to the 
extreme value is still open. 

We know that the frequencies and damping times of the quasinormal modes are entirely 
fixed by the black hole, and are independent of the initial perturbation. It has been shown 



that there is a second order phase transition in the extreme limit of black holes [28 , p9| , 30] . 
This result has been further supported in a recent study for the charged AdS black hole 



31| , p2[ . Whether different properties of imaginary frequencies with the increase of the 
charge Q reflect different phase characteristics is a question to be answered. We are trying 
to improve our numerical code to increase convergence rate and precision. These results 



can be compared to those of |33| where it has been shown that in the Kerr solution, when 
the rotation parameter a of the black hole approaches the critical value, the quasinormal 
modes decay time gets longer. These results support the claims of the present paper. We 
will address this problem in detail in a forthcoming paper. 



4 Conclusions and discussions 



We have studied the wave evolution in Reissner-Nordstrom Anti-de Sitter spacetimes and 
revealed consistent but interesting results. The radiative tails associated with a massless 
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scalar field propagation on a fixed background of an AdS black hole experience an oscillatory 
exponential decay at the black hole event horizon. This result of exponential decay supports 
p0| and is similar to that on the black hole horizon in de Sitter case because of the similar 
behavior of the exponentially falling off of the potential at event horizon. When one considers 
the conclusions of Ching et al flB] , this late-time tail is not surprising. Accompanied with 
the exponentially decay, there is a special oscillation in the tail in linear analysis in AdS. 
This can be attributed to waves bouncing off the potential at large r, because the potential 
diverges at infinity in AdS space. The property of the late-time behavior on the event 
horizon of AdS black hole can help us better understanding the stability of Cauchy horizon 
in AdS spacetimes. 

We have also learnt an object-lesson of the quasinormal modes on the background of 
RN AdS black hole. For small values of Q the picture we obtained is consistent with that 
derived in J22|. The larger the black hole charge is, the quicker is the approach to thermal 
equilibrium. The consistent picture of the quasinormal modes depending on the multipole 
index I has also been illustrated. Increasing /, we obtain the effect of increasing the damping 
time scale and decreasing of the oscillation time scale. 

When the charge in background RN AdS black hole is large enough, we get an opposite 
result for the characteristic of imaginary frequencies. This "wiggle" of the imaginary frequen- 
cies as the charge increases on the background agrees with earlier results of normal-mode 
frequency study for RN black hole. Since the type of quasinormal mode is determined by the 
background spacetimes and it has been found that a second-order phase transition appears 
when the RN AdS black hole becomes extreme, we speculate that the different behavior of 
the quasinormal frequencies may reflect the characteristic of two different phases. Further 
study on this subject is called for. 
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A Convergence of the Code 



As shown in expression fll7j) , the local truncation error for ^jy scales as e 4 , assuming that 



the exact function has a Taylor expansion in the null rectangle. Although we have little 
control over the global error at u = u. max , we can make a rough estimative. The first point 
in the line u = u max , which is ty(u max , vi), is obtained using only points in the line v — V\. 
Since there are 0{e~ 1 ) points in this line, the global error in ^i(u max , vi) should scale as e 3 . 
On the other hand, since there are 0(e~ 2 ) grid points, the last point in the line u = u max , 
which depends of all the grid points, should scale as e 2 . 

We can improve on this estimative calculating convergence curves. In these curves, we 
plot the value of at a given point as the grid spacing e is reduced. In the figures 8 and 9 
we show some curves for three value of the charges in two points. The parameter N is the 
number of points in the line u = uq. It is therefore inversely proportional to e. 
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Figure 8: Graph of for v = 2.9824, r + = 0.4, / = and several values of Q. 
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Figure 9: Graph of for v = 9.6157, r + = 0.4, / = and several values of Q. 



These diagrams clearly show numerical convergence, and a quite fast one. This indicates 
that the global error indeed tends to zero as we decrease the grid spacing, even in the limit 
when Q gets near the extreme limit. 
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